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Statistical Properties of the Time Histories of Cosmic Gamma-Ray Bursts 
Detected by the BATSE Experiment of the Compton Gamma-Ray Observatory 


Final Report 


The main scientific objectives of the project were the following activities: 

(1) Calculation of average time history for different subsets of BATSE gamma-ray bursts. 

(2) Comparison of averaged parameters and averaged time history for different BATSE 
GRBs sets (by Peak Flux, Hardness Ratio, etc.). 


(3) Comparison of results obtained with BATSE data with those obtained with APEX 
experiment at PHOBOS mission. 

(4) Use the results of (l)-(3) t0 compare current models of gamma-ray bursts sources. 

To obtain these results, special statistical studies of 338 cosmic events of the Second BATSE 
Catalog have been done. The corresponding conclusions follows: 


(1) The new method was developed by averaging time profiles of separate events when they are 
aligned around their principal peaks. This method was compared with other methods developed 
for similar analysis. An agreement was found in similar datasets used in these tests. 


(2) Using the average time histories at different energy channels on 1024 ms time scale, the following 
features were found: 


- The average time history of GRBs is non-symmetric: from the level of 0.25 of its maximum 
height, the rise front of the average curve of emissivity (ACE) lasts about 3 sec and back slop 
about 5 sec. The average Hardness Ratio has a broad maximum along the rise front of the 
ACE, then decreases along the back slope. 

- From the comparison of ACE obtained at different channels, we find that the averaged 
peaks of GRBs are broader at soft energies then at hard energies (softness/duration effect). 
Using the peak shape of discriminator No. 3 as the standard, we find the best stretch 
factors K to fit the peaks in discriminators No.l and No.2 , as the following low: 
K(E) = (173 (keV}/E) 0 23 . 

- We find no evidence in this analysis for time-dilation in the ACE of weak events relative 
to bright events, when those sets of events are selected by the photon peak fluxes, as 
presented in the Second Catalog. The ACE at energies 50-300 keV for dataset of 143 bright 
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events with F>1 ph/cm 2 s and 179 weak events with F<1 ph/cm 2 s correspond to similar 

curves. In the more conservative case with the gap between peak fluxes of bright and weak 

GRBs, the average curves have been shown to be similar as well. 

— We find the difference in average curves of Hardness Ratio for the sets of bright and dim 
events selected by the photon peak fluxes. Bright events have larger Hardness Ratio then 
weak ones not only averaged over time but at each point in the averaged time profile (effect 
of the brightness/hardness correlation). 

(3) The comparison of results obtained with the BATSE data have been done with those obtained 
with the APEX/PHOBOS data. The perfect agreement is found between them. 

(4) The following conclusions have been obtained for the theoretical models of GRBs. 

— Similarity between all outbursts and of all outbursting objects was found. This fact rules out 
all galactic models, which attribute GRBs to different populations of neutron stars. Another 
modem galactic model could also be rejected, which attributes all GRBs with the single class 
of neutron stars, but which are thought to emit two different types of strong and weak 
outbursts. On the other hand, the perfect similarity between average curves for bright and dim 
events makes the cosmological paradigm of GRBs quite unprobable also. 

- Average evolution of outbursting objects during a burst shows non-symmetrical rises and 
decay fronts on the average, the developing phases of GRBs before the maximum are shorter 
and harder then decaying phases after the maximum. The picture of outbursts is not 
time-reversible. Some sort of physical evolution of the outbursting object does occur during 
a burst, which leads to these features. 

— Theoretical models of GRBs agree with the brightness/hardness correlation of 
GRBs in two cases. In the first case, the model might explain the difference in brightness 
from different distances of sources, assuming that some evolution takes place. There is 
another possibility also, when the all range of brightness of GRBs is associated with the 
range of intrinsic luminosity of sources, which are at the same distance from the observer. 
Further study is necessary to resolve these possibilities, but all models with the standard 
candles could be ruled out due to brightness/hardness correlation of GRBs. 

The main results of the project were reported in: 

- 2nd Huntsville Symposium on Gamma Ray Bursts; 

- 30th Scientific Assembly of COSPAR; 

- 17th Texas Symposium on Relativistic Astrophysics. 

The main results were presented in the following publications: 

[1] Mitrofanov et al. 1994, In Proc. of "Gamma-Ray Bursts", AIP 307, p. 187. 

[2] Mitrofanov et al. 1995, subm. to Astrophys. J. 

[3] Mitrofanov et al. 1994, In Proc. of 30th Scientific Assembly of COSPAR, to be published. 
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ABSTRACT 

The results of new method of statistical studying of APEXJPHOBOS and BATSE/COMPTON cosmic 
gamma-ray bursts (GRBs) are presealed, which are baaed on the averaging of time profiles or separate 
events. Using time histories at different energy channels, the general features are found of averaged 
fUix/spectnini evolution of GRBs. The comparison is done between fWx/spectnrm evolution of bright and 
dim events: while no cflffereace* were found between averaged curves of fluxes variability, the comparison 
of avenged evolution of their hardness ratios pointed out the effect of haidncss/brightness correlation. The 
consequences of the statistical results are discussed for two principal alternatives: for the cosmological 
model of GRBs and for the models which associate them with the Galaxy. 

INTRODUCTION 

GRBs are known to have very random time histories [1], Even assuming that sources are standard 
candles', one should suppose that some randomizing factors ire acting during emission, providing very 
large changing between different events. The comparison of Inch vidua] events hardly might resolve such 
systematic effects, as differences between averaged flux/spectnun variability. 

On the other hand, the averaged flux/spectruin variably has the physical sense, which represents the 
general properties of omb orating sources and basic features of bursts mechanisms. The future theory of 
GRBs, as the single phenomenon, has to explain the averaged flux/spectrum variability. Moreover, ai the 
present Level of knowledge of GRBs* some differences might be expected in the flux/spcctrum variability of 
bright and dim events. There are lot of new or revteUed models of GRBs, which deal with some of these 
effects; e g. the model of two types of outbursts from close galactic neutron stars [2], or the model with 
two ootborsting populations [31. or the cosmological model with necessary predictions of time stretch and 
frequency shift |4), 

Those general observational features hardly can be found in the study of particular GRBs, but they might 
be found In the special averaging of the data available. And* provided some of them would be found, (he 



origlu of ORBs couW become dear, because particular effects are associated with specific models 
Thexeiore» special statistical analysts of the total set of available ORBs is now very actual. 

Together with new models, new methods of statistical analysis have been proposed to resolve the ongm of 
ORDs, In particular, the methods of averaging of time histories have been developed 15-9). Ajxrther known 
inethreto wore baaed on the wavelet analysts of GRBs. which provide the averaged activity curve* for 
different acts of selected evema ( 10 - 12). The separation of GRBs into two different classes of variability 
has been proposed (13|. This studying Interferes with the method of estimation of time duration of GRR« 
(14), which permits to resolve two possible classes of GRBs with different averaged duration, namely 
smaller and larger of 2 seconds [15). One might conclude that some “new statistics' of GRBs arc appearing 
now, which represent the general features of flux/spectnim variability. 

This paper presents the results of the particular method of new statistics' of GRBs, which Is based on the 
averaging of events aligned around their principal peaks. The studying focuses on the general features ol 
ffuK/spectrutn variability and gives the way of simpfc comparison between bright and dim eveiusi. 

THE AVERAGING OF FLUX/SPECTRUM VARIABILITY OF GRBS 

To resolve real general features of flui/spectnim variability of different GRBs, one needs to develop 
method, which excludes a$ much as possible the influence of their different brightness and their individual 
time histories. The methods have been developed [5-9), which respond to this demands. They include the 
following procedures: 

1) For averaging of an selected GRBs, a spectral range should be selected, where their time histories will 
be averaged. It is possible, that different spectral ranges could he used for alignment of time histories, and 
for averaging of aligned events. 

2) A unique time scale should be selected with time bins, which ensures acceptable signal to noise ratio at 
principal peaks of all time histories, even for dimmest events. Bin* at selected time scale have to 
correspond to the actual exposures of each events in the selected energy range. 


3) Along time history of each hurst, rise front, back dope and the unique time bin of primary peak 
should be selected. At the primary peak the maximal number of counts should be determined After 
that, normalization of time history should be done, when actual counts C\ at all time bins t* are replaced by 
dimensionless values t x t=Cyc: mw , 



flux/specirum variability of GRBs: 
(I) and (II) - Averaging I and II. 


4) The first method of averaging of time histories 
{Averaging f) provides ihe Joint statistical study of time 
intervals wtih increasing and decreasing counting rate at 
time intervals of particular events, which have Che similar 
dimensionless fluxes. For selected intervals of relative 
flux F, the averaged durations tf^p 35(F) of all lime bins 
arc obtained with corresponding fluxes for 
rise fronts (RF) and back slopes (BS). separately, for 
all contributing GRBs (Pig. 1), the values ortfcp(Fj) 
and tgs(Fi) represent the averaged curve of flux 
variations for selected set of GRBs. The errors of 
estimations of those time intervals represents 
the errors of used relative fluxes of particular events. 


5) For another procedure of the averaging, all normalized time histories should be aligned around the same 
time bin, placing their peaks at the single bin to- For time bins tj Of common time scale before and 
alter Iq, Die values of normalized fluxes fj are averaged over the selected set of bursts. The second 
procedure of averaging G Averaging it) builds the dimeittlocilesa value F] along Ihe Ume scale, equals 10 ti*c 
sum of f j of each contributing events, efivided by their number for this time bin (Fig. 1). The values F| 
along time scale represent the averaged curve of cnUssivity ( ACE ) for selected set of GRBs. The crior ol 




each value Fj is associated with statistical taota of measured numbers of counts Cj. Theca errors r&proecnt 
the uncertainty of ACE, which la resulted from the uncertainties of particular inaaBuremcius. which is 
much smaller for bright bursts then for dim events, 

6) Estimation of the hardness ratio (HR) of detected counts Is the moat simple way to Study the evolution 
of energy spectra along time histories of CRB*. For the Averaging J. the moan values of hardness ratio 
could bo e st i mated Tor each separate lime intervals Irj^bs^i) The >' re P r8 ““ u 0,0 av * ra S e<{ HR f or 
selected Interval of relative fluxes at rise front* and back slopes, respectively. The averaged ciuve of HR 
(AC HR) could be created along (he time scale of the Averaging //. as the ratio of corresponding hard ami 
soft ACE. averaged around aligned peaks at the integral spectral range. 

AVERAGED FLUX/SPECTRUM VARIABILITY OF GRBS 


■Study of APE X/PHOBQSriala 

The Russian-French experiment APEX on the interplanetary spacecraft PHOBOS-2 detected 62 GRBs 
during 8 months of flight from July 1988 till March 1889. The tostnimcm used a 10 x 10 cm CsI(Tt) 
scintillator to measure time profiles In the broad energy interval 120-1420 keV and energy spectra in 108 
channels from 65 to 9000 keV. The totality of 48 APEX/GRBs with good {demented data have been 
loaded to the APEX data base for further statistical study 15-71. For the averaging, toe APEX/GRBs are 
divided into subsets corresponding to the 2 sub-sets of 19 bright and 29 dim events with V/V max at sub 
ranges of (0.0 < V/V^ < 0.3) and (0.3 < V/V^ < 1 .0), respectively. 

The Tbble 1 (6] represents the comparison between toe mean durations Irf ( bS fi* ^ interval of relau ve 
fluxes [0.35 - 1 .01 at rise fronts and back slopes, respectively (1/2 of the duration of peak intervals f mux = 1 
1/2 are included into RF and BS). For relative fluxes greater then 0.35 one codd be sure that all me events 
considered, the total number of counts contained in the broad channel 300-1000 keV would exceed the 
level of 3 sigma od background fluctuation*. The Table t ilao give* the averaged v*lus of HR, defined as 
HR*F(300-1000 keVyF(100-300 keV) corresponding to the Intervals of relative fluxes 10.35 0 70JaiRP 
and BS. To estimate HR. summation has been done of all contributing particular energy spectra, which 
have been normalized to ensure the equal contributions of all selected events. 

TAB! F. 1 Averaging I [61 for intervals of similar relative fluxes at rise fronts and back slopes 
of APEX/FHOBOS GRBs. 


Selected set of GRBs 


Rise front, (RF) 


Back slope. (BS) 


Set of bright events 
0.0< V/V^ <0.3 


Set of dim events 
0.3 < V/Vmaj, < 1.0 


tRp (0.35-1.0) = 6.2±0.8s 
HRrf<0.35-O.70)=0.40±0.OJ 
IrP (0.35-1.0) - 3.1±0.9s 
HRrf(0.35-0.70)=0.25±0.03 


l BS (0.35-1.0) = 3.6±0.9s 
HRbs( 0.35-0.70)=0.3 l±0.01 
tp S (0.35-1.0) - 3.6±10s 
HR BS (0.35-0.70)=0 24±0.2 


Rom the Table 1 one might conclude the difference between avenged RFs and BSs of bright GRBs and 
their symmetry For dim events, and the difference between tpj; and HRjjp of bright and dim events. 


Another studying of fhix/spcctrum variability of APEX/PHOBOS GRBs using (he procedure of du- 
Averaging // permits to see the general properties of toe evolution along the bursts. The averaged curves 
of tlux/spectnim evolution tor bright and dim sub-sets of GRBs are presented at the Figure 2. Both ACE 
curves have one peak shape with smooth back slopes. Those back slopes are very similar for bright and 
dim events. Applied statistical criteria do not show any significant differences between them [5.6). 'fhc 
difference at the rising parts of ACEs, with curve for bright tub-set going higher then one for dim events, 
is probably caused by toe selection effect: For weak events the triggering happens earlier for sleep rise 
fronts then for gently sloping ones, therefore dim events might be in a deficit before toe maximum. 
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Figure 2. Averaging II [7J: ACE (thin line®) and ACHR (thick Unas) for bright sot 

0.0< V/Vm «<0. 3 (aolid Unc, 1 9 events) and dim set 0.3<V/Vtna*<1.0 (dotted line, 29 events) for 

APEX/PHOBOS GRBs.. 


On the other hand, after the maximum, there is no influence of the triggering at all, and shapes of ACE 
curves represent the real cose®. At the TaWo 2 the results of the Averaging U ore presented for sets of 
bright and dim APEX/GRBs, which compare the averaged fluxes and HR along the ACE curve. 

TAB 1 JE 2 Averaging II [7] along lime curves of APEX/PHOBOS GRBs 

Selected set of GRBs 1 * before t roax Bln at 1 s after t max 

Set of bright events F=0.53±0.01 F=1.0 by definition F=0.53±0.02 

0-0<VA/ max <0.3 


Dim events with 

o.3 < vr/ nux <i.o 


HR=0,37±0,02 HR=0.45±0.01 HR=*0.29±0.02 

FM).49±0.03 P-1.0 by definition F=0.46±0.02 


HR=0.27±0.03 HR«0.25±0.03 HR=0,30±002 


While no difference Is seen between ACE curves for bright and dim events, the hardness at fronts and 
peaks of bright GRBs is much higher. However, the APEX/PHOBOS experiment have provided very 
limited sets of GRBs, and much better statistics of BATSB/COMPTON Is necessary for statistical 
studying of averaged flux/spectnun variability with higher significance. 

Study of BATSE/CQMPTON data 

For statistical analysis of fluxApectnim variability of B ATSE/COMPTON GRBs, 260 events have been 
used [8,9) at the First BATSE Catalog [14]. The events have been excluded with short lime duration 
* 90^1 see, and other ones with missed data; the total number of excluded events was 52. Die used set of 
208 BATSE GRBs corresponds to data recorded from April 19, 1991 till March 5, 1992. For each event 
of 260 GRBs from [14], the data rile has been used of measurements of Large Area Detectors (LAD) with 
time resolution of 1.024 s in four discriminator channels No.l( 25-50 keV), No.2(50-100 keV), No 3(100- 
300 kcV) and No.4(>3(X) keV). To create the time history of an evert, data has been taken from LAD, 
which had the largest counting rate. From the time profile of die continues measurements, tltose time bins 
have been selected fear time history of an event, which correspond to Intervals of ±tnax(t 9 Q, 20 s) around 
the moment of burst triggering Measurements outride those intervals have been used to interpolate 
and subtract background during an event 
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The Averuing I study of the B ATSE/COMPTON GRBs corresponds to much better statistics, then in the 
case of the APEX data. The first results of this analysis [8,91 are presented below. 



The averaged curve of emitstvUy AC&>3) curve for the frit set of BATSE GRBs |14) is shown at 
figure 3 for discriminator channels No. (2* 3) at lie energy i**e or 50-300 keV. The curve has oo^peak 
structure with non-syrmnetrical from (F) and decay <D), wMch have widths above (he level of 0.1 equal to 

t fd = ( i + f i (F,D) ) • 1024 * (u 

forf^ 0 ^ 

Tlie total width of peak above the level of 0.1 equal* to Tp- Tp + Tjj. Corresponding values of widths we 
presented ill the Table 3. It is evident from this data that averaged front of ACE is much shorter then 
decay. TO itlusinic non-symmetry of the peak, the lime- reversal decay la shown at die front side (Fig 3). 

TABLE 3. Averaging 11 [9]: Widths of fronts, decays and total peak of ABC above the level of 0. 1 
for full set of BATSE GRBs |14] (sum of discriminators No.(2+3) at energies 50-300 keV) 


Selected set of GRBs 

Tp (s) 

T D (s) 

T P (s) 

Full set Of 208 events 

2.36*0.02 

3.83*0.02 

6.19*0.03 

81 bright events 
with flux >lph/cm 2 s 

2.65*0.01 

3.82*0.01 

6.47*0.01 

127 dim events 
with flux clph/cm^s 

2.11*0.02 

3.78*0.03 

5.89+0.04 



Figure 3. The ACE(2+3) curve for full set of B ATSB GRBs 1 14| for discriminates No. (2+3) at 
range 50-300 keV [8,9]. The Ume-revcrsal decay of the main peak Is shown with dash line. 

Time histories at discriminations No.l-No.3 could be used separately to build corresponding ACE's The 
estimations of widths of fronts, decays and total widths of ACE's peaks at discriminators No. 1-No. 3 arc 
presented in the Table 4. From the comparison of ACE curves obtained at different channels, one might 
conclude that averaged peaks of GRBs are broader at soft energies then at hard energies. Using the peak 
shape at the discriminator No.l , as the fitting model, one might find IIS best stretches to Ot two other peaks 
at discriminator No.2 and discriminator No.3. The corresponding stretch coefficients K(2.3) are presented 
in Table 4. Using those values of K, one might Interpolate the stretch/ energy dependence within the range 
from 25 up to 300 keV, as the following: 


K(E )„ ( 170*wV/ e )« 


with the best fitting values of parameter a equals to 0.24*0.03 for reduced c/iA=0.5. 


( 2 ) 



1 * | ARr f. 4 Averaging U [9J: Widths of fronts, decays and total peaks of AEC above thej level of 

0.1 for discriminators No.l, No.2 and No.3 for fall set of BATSE GRBs (14] i 



Discr. No. 1 

Discr. No.2 

Discr. No.3 

Tr(s> 

3.07*0.03 

2.52*0.02 

2.02*0.02 

Td(s) 

4.76*0.04 

4.01*0.02 

3.10*0.02 

Tp(s) 

7.83±0.05 

6.53*0.03 

5.12*0.03 

K 

1.43*0.08 

1.27*0.05 

1.0 by defini tion 


The ACE curves for BATSE GRBs could be compared with similar avenged curves previously obtained 
for GRBs delected by APEX experiment. Data of 48 APEX GRBs at the energy range off 1C 0-1000 1 keV 
could be compared with ACE for 205 BATSE events far the discriminator No. (2+3) (FI 5. 3). Good 
agreement between averaged curves of APEX and BATSE events shows that they represent ieai physical 
evolution of emissivity GRBs and faintly depend 00 type of instruments used. 

The fact that AEC depend on the energy means that some kind erf spectrum evolution shoul 1 take place 
al ong me averaged light curve. This evolution could be resolved from the comparison of 1 1 e averaged 
values off Hardness Ratio (HR) for different bins along ACE. To estimate averaged HR, ( Iscriminator 
channels No.(2+3) could be used to align the primary peaks, and averaging should be done tased on this 
alignment of time histories at channels No.2 and No.3. The ratio of ccnespowjbig curves fer #3 and #2 
provides the averaged curve of hardness ratio variability (ACHR) (Fig. 4). 

The time profile of ACHR shows the broad maximum at the hoot of ACE curve dialog about 0 seconds 
before the main peak. After the peak, the averaged hardness ratio decreases with the dccreasinf averaged 
flux along the decay. The broad mTimum of hardness ratio at the rise front of averaged peak < if ACE 

means that usually peaks at higher energies take place before peaks at medium and low energies. This 

conclusion, based on the present statistical analysis of BATSE GRBs, confirms the pervious o nervations 
of many particular GRBs, when peaks of hardness ratio have been found at fronts of peaks [ 1< ]. 



Figure 4. The averaged curve of hardness ratio (ACHR) (8,9) for full set of BATSE ( >RBs [14], 
with hardness ratio defined as (100-300 keV)/(50-100 keV) (8.9). 

To compare bright and dim BATSE GRBs, two sub-sets have been fanned from the I'irst Catalog 
database. As hright and dim GRBs, those sub-sets of events have been selected, which have 1 > 1 ph/cm s 
and F<1 ph/cm 2 s. respectively. Taking into account that the time scale far averaging has resolution of 
1.024 s. only those GRBs have been used for comparison, which have duration t90>l s. Ttyre were two 
sub-sets of 81 bright and 127 dim BATSE GRBs, which have been used to compare ine averageo 
flux/spectrum variability. | 




For the sets of blight and dim GRBs, the ACEa curves for discriminator channel* No, (2+3) arc presented 
in Figure 5. Both curves have one-peak profiles, and no significant difference is evident between them. IIjc 
width of rise fronts of ACB peak of bright events Is slightly longer then Uie similar width of dim events, Dul 
both back slopes have very similar duration (Table 3). 



Figure 5. ACE curves [8,9] for sets of bright and dim BATSE GRBs [ 14] for discrinunaiors 
No.(2+3) at enery range 50-300 keV. 

The ACE curves of bright and dim sub-sets could be compared more accurately by the using of central 
peak for ACE of bright sub-set, as a stretching fitting model for central peak to ACE of dim events. 
Taking into account only those time bins of main peaks, where ACE curve lies above Uic level of F> 0.1, 
foe best fit of ACE^) corresponds to the oanwing factor K*0.97 ± 0,05 of ACB< brt ^0 with reduced 
cftA-2,80. 

Thus, from the comparison of bright and dhn sub-sets of BATSE GRBs, the conclusion should be drawn 
that there is no statistically significant tflffiercnce between them in the averaged cmissivity curves. All 
proposed methods of comparison do not show that ACE for dim set is broader liven ACE lor bright set. 
Moreover, there is even some evidence that It Is narrower then ACE for bright events. 

Another comparison could be done between sub-sets of bright and dun BATSE GRBs, which is associated 
with their mean hardness ratios. Two ACHR curves were calculated for two selected sub-seta ol bright and 
dim events (Figure 6). Those two curves have similar shapes wilh broad maxima at fronts and peaks 



hardness ratio defined as (100-300 keV)/(50-l00 keV) [8,9]. 

of ACE. The most important finding of this comparison is Ihe difference of die levels of those two curves; 
along the full all time scale the curve ACHrO**^) goes well above the curve ACHRW* m X The 
difference between AC HRs curves for bright and dim events could be interpreted, as the evidence for the 
effect of brightoess/hardness correlation of BATSE GRBs. 





CONCLUSIONS 


The main results of the statistical studying of fluxArpcrtnim variability arc the following: 

i) Bright events have longer and harder averaged rise fronts then back slopes. Dim even la have 
more symmetrical averaged rtee fronts and back slopes, both in the averaged duration and in the 
hardness ratio (Averaging lot APEX/PHOBOS GRBs [5,6]), 

11) Averaged rise fronts are longer and harder for bright events than for dim events. Averaged back 
dopes of bright and dim events have similar time duration and hardness ratio (Averaging / of 
APEX/PHOBOS ORB* [5,6]). 

ill) Averaged curves of ccnissivity for both bright and dim events are very similar, especially 
at their decays. Averaged curve for HR of bright events goes higher then for dim events [Averaging 11 
of APEX/PHOBOS GRBs [7J and BATSE/COMFTON GRBs [8,9J). 

The first evidence of the effect of brightness/hardness correlation (see (ti) and (lii)) luts been found using 
the APEX experiment data, when two averaged energy distributions have been compared, which 
summarized particular energy spectra measured at peaks of strong and weak detected events [17]. Later, 
when the method of averaging with peaks alignment has been developed, the presence of 
brightnessfoardness correlation has been demonstrated for APEX events more dearly [5-7]. Independently, 
the effect of brightness/hardness correlation has been resolved la the limited set of the first 126 BATSE 
GRBs [1 81. Finally, the studying of the set of 205 GRBs from the First Catalog of B ATSE had shown Utc 
presence of hardness/brighiness correlation [8,9]. Using another set of bright and dim BATSE GRBs, 
larger spectral hardness was found for bright events in comparison with dim bursts [19], 

CosmctogY Model 

Developed method of statistical studying of time/energy evolution of GRBs permits to compare bright and 
dim events according to the predictions of cosmological model. To explain observed value of averaged 
<yrv wa > for the first 1 18 BATSE events, the cosmological model has been developed [2D], which works 
with sources of GRBs, as standard candles. For the best fit of the observational data, die model 
corresponds to a red shift about unity for a most distant objects, but the factor of 2 is possible for other 
marginally acceptable fits. 

Assuming the validity of this model, the following testing procedure could be used based on Uic metixxJ 
proposed. Tlie sub-set of brightest events is selected, which corresponds to the most close sources with die 
smallest rod-shift, one might take those 38 events, which have 0<V/V, niU <0.05. According to the 
cosmological model (20], their mean red-shifts could be estimated, as £fcr*0 ) t. 0.29 and 0.50 for low z 
fit, best z fit and high z fit, respectively. On the other hand, to select the most distant observable sources 
with the largest red-shift, one might use the comparable number of 42 dim events with V/V max from die 
interval 0.CxV/V nux <L0. According to the cosmological model, for those cfcm events one estimates 
z t jj jn =l.24, 1 36 and 1.70 for low z fit, best z fit and high z fit, respectively. 

According to the cosmological model [20] and the width/energy dependence (2) (sec Table 4). the 
theoretical evaluations have been done of the total stretch factor, which provides the best fit of ACE cd 
dimmest events at discriminator No. 2 by the ACE for brightest events at discriminator No. 3. 'Jhosc factors 
were estimated, as 1.42 for low z fit, 1.53 for best z fit and 1 .83 for high z fit [9]. On the other hand, the 
best fitting value of the stretch factor [9] equals to 1.04*0.1, which is far away from the theoretical 
predictions. So, the averaged ABC curve for brightest events does not differ so significantly from the AEG 
curve for dimmest events [8,9], as it could be expected according to the cosmological model [20]. 

Galactic Models 

The averaged emlssivity curves AEC are very similar for sub sets of bright and dim events, This fact does 
not agree with models, which attribute GRBs to different populations of galactic neutron stars (e.g. [3, 
21 1). These classes arc thought to belong to the close vicinity of galactic disk and to the extended galactic 
halo. Gore sources might be responsible for dim events only, because more distant outbursts from iliusc 



population, winch would dominate at low galactic latitude*, are not observable. On the contrary, distant 
neutron stars in the halo might be attributed to bright CRB». Talking Into account the ratio at bright and 
dhn events -10 2 - 10 3 . and Uic difference in their distance* - JO 3 , one might deduce the difference o( 
intent at radiated energy by Otoac two classes - 10® - 10®. Such a large difference in the radiated energy 
pot ota out. that automating activities of those two classes have very different physical nature 

Another modem galactic model attributes all GRBa with the stogie class of neutron stars in galactic disk, 
but which arc thought to emit two dUXerent types of strong and weak outbursts [21. Close members of this 
population are thought to be visible either through events of the strong class, observed as bright C-RBs, w 
(trough events of weak class, observed as dim GRBs. More distant pan of the population is drought to be 
detectable through events of the strong class only, observed as medium GRBs. Therefore, ior Uic model 
[2], performed statistical comparison of bright and dhn events corresponds to the comparison of iwo 
different classes of outbursts from the close stars, i.e. classes of strong and weak events. 

Therefore, taking into account the similarity of ACE's for bright and dim sets, the model with different 
populations of GRBs sources does not seem to be very probable, as well as the model of single population 
of sources with different classes of outbursts. 

While the perfect similarity has been found between the averaged emissivity curves of bright and dim sub- 
sets of GRBs, as the evidence for the homogeneity of radiating objects, the difference in Uic averaged 
hardnesses has also been seen between those two sub- sets, which points out on some intrinsic difference 
between emitters of bright and dim events. 

The model of extended galactic halo population (e.g. see [221) agrees with the hrightncss/haidness 
correlation of GRBs in two cases. In the first case, the model explain* the difference in brightness by 
Afferent distances of source* within the halo, assuming some evolution of non-standard sources “Jong the 
run away path from the center region to (he periphery (e,g. see [23. 24». For those models, the mam 
differe nce between brighuras of GRBs in -10* times is attributed mainly to the broad spread of their 
distance* 3-300 kpc. The change of sources could be resulted Aon some sort of evolution of outgoing 
neutron stars during -300 (v^lOOO lcm/s) Myra. In this case the effect of brtghlocss/hmdness correlation 
might correspond to some physical correlation between inverse distance and hardness of gamma-rays. 


There Is another way to implement the idea of non-standard sources into the model of extended galactic 
halo, when the all range of brightness of GRBs is associated with the range of intrinsic hinrinosUy at 
sources about ~ 10* times, which all are atthe comparable distances -30-300 kpc from the observer. The 
reason is unknown for the absence of GRBs from the inner regions ofihe halo, and neutron stars have to 
preserve the ability to radiate GRBs during all maximal age of -10* Myix In tins case the effect of 
brightness/haidiiess correlation might be associated with some intrinsic correlation of radiated power with 
hardness. 
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ABSTRACT 

Find rosuRg urc presented fi x*m the analysis of 200 HATSE gammi^ray 
bursts (GltUa) 1 using & meUiod l,$ whereby all event* ore synchronized and 
averaged around the bine which arc the brightest time interval* of each of them. 
Fbr the averaged time history* a difference i* found between the rise front and the 
back slope, and good evidence is found for the presence of hard-tosoft spectral 
evolution. We compare sub-seta of “strong" and “weak” events and find no 
evidence for time dilation in weak GRB», as would be expected in coemologica] 
GRB models. On the other hand, for the strong events the averaged hardness 
ratio is found to be larger than for the weak GRB*. 

1. INTRODUCTION 

A* discussed in many contributions to this Workshop, our present knowl- 
edge of GRBs is full of apparent contradiction*. One of these in the contradiction 
between the perfect isotropy of GlLlfe on the sky 4 **, which corresponds to a very 
itomogcncousi angular distribution, and the value of {V/V mhX ) G.32 4,8 , which 
indicates that GRBs have a non-homqgeneous radia] distribution in Euclidean 
apace. 1b resolve this contradiction one might suppose, for Instance, that GRBs 
axe generated in the Extended Galactic Corona 4 , Or that they belong to a cosmo 
logic*] population with redeidft * 1 T . A definite decision between cosmological 

and nou-cosmological model* would he a crucial step in understanding the origin 
of GRBa 

The method of averaging time histories cif cosmic GRB* around the bright- 
est time interval was devdoped for analysis of event* detected in the APEX 
meperiment of the PHOUOS Mission 2 . A set of 48 event* ww divided into two 
subset* of ‘’strong" and “weak" events* and their averaged rima histories vrere 
directly compared. While different* were found In the rising parts, the hack 
blopes were identical within statistics. The difference* in the riling parte could 
result from a systematic difference in triggering moments for strong and weuk 
GRBs. On the other hand, tlie similarity of l>ack dopes indicated that there is 
no cuwndogic*il stretch of weak events with respect to strong GRBs. 


1 also University of Alabama in Huntsville 
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Jl. AVERAGE TIME HISTORY AND HARDNESS RATIO 


The data from the BATSE instrument on the Compton Gamma Ray Ob- 
wjrvatory provide an excellent opportunity to study the averaged time histories 
of GRBtt, in particular to compare subnets of strong and weak events, with much 
better statistics. Wo prauiat hsrsin prcUmiiiaiy rwwilu of this comparison. 

A total of 260 GRBs from the first BATSE Catalog* have been studied so 
tar. For each event we used ths DISCLA count rate time tmlorim (four broad 
energy channels with a time resolution of 1 1 ) from the large area detector which 
hsd the highest peak count rate. Those are independent of Uni instrument 
trigger* so that one has (dr every event a continuous tone history containing 
both pro- hurst and post -burst background W® defined each burst interval as 
extending from Ti — AT to To + Af, where Ho is tlie burst trigger time and 
AT *• uiajc(Tao* 30 *), where T§* is a measure of burst duration 1 . Background 
during thu burst interval wvut interpolated using the measurement* before *rtH 
after this interval, and tho difference between The measured actual counts and 
tiie interpolated background curve constituted the burst time history. 

rhe results we report arc derived from two brood discriminator channels 
#2 (50-100 keV) sad #3 ^100-300 keV). Fbr each GftB time history we deter- 
mined the brightest time bin where the maximal number of count* w« npmrtM 
in the sum or channels #2 and The time histories (with statistical uncer- 
tainties) Ln the separate channels were then normalized to the count rate in the 
brightest bin. The result at a standardised set of GRB time histories, with time 
resolution of 1.1)24 a and maximum magnitude equal to 1. 

With this data typo, all errant* with duration *hnrirr than 1 a are unre- 
solved, their time historic* falling within one or two bins, depending on their 
pLid*c with respect to the instrument dock. At tfosstago, ail events with 1^* < 1 
:i were therefore uxcludod from ths analyst*, fu the imafc analysis ptiase* we plan 
to iwc data with a time resolution of 64 ms, allowing use of the shorter event.* 
Tlie averaged time history in the sum of discriminator dtuuuieU #2 and 
#3 (60-300 heVJ i» pretented in Figure 1 for 206 GRBs with T*j > 1 $ from tl»e 
hist 15A1SE eiUAkHt*. Betdd w Uinshutt Vinau, a vUm f««»m kW. CTnSftlng 

<XKilii not be iwcJ Juu to data gap. The wwacuoo oi Umkw svmuU *1*4 UM WW All* 
the final reflfdt: the sot of 206 GuBu fully represent the continuous observation 

S criod of BATSE from Ifi April 1991 to 5 March 1992*. Tl>e time scale of ail 
gures corresponds to numbers of 1.024 * bins before and aRer the brightest 
one, which vc defined uo wero. Within the time interval from bin -lii to bin 20, 
the total number of contributing tfvoaU (i.e., with significant count role above 
background) is more that 90% of the total number in the averaged set. The 
mean dura clou of ibe avuiugnd iuus p*«jj(Uo a* the lcvri of 0.36 «f ito inuumuxn 
height is about 7-9 a, which corresponds to 2-3 s of rise (rant and 4-^5 s of back 
dope (with uncertainties of about 1 »)« The difference between the averaged rise 
(rout and b**4t alope is dluntratod by the dashed curve In Fig. 1, which shows 
the time- reversal of the averaged light curve. 

GRBa are known to have broader piyJc* of time histories at low energies 
relative to high taiagic*. La Fig. 2 vc compnro tlm art) rafted time bwinry r>f 
the full sot of BATSE GRBs in discriminator channel #2- aligned around the 
Uiigtowt VJa { ~ —UV. l:»rrmM<itnr« in channel #3 aliened 

unwind tho brightest time bin in channel ^3. It is easily seen riwil the softer 
uuurgy channel indeed has a broader peak in it* averaged tiaie hmtory. 

Tho Averaged tum. history of BATSE CRD* may be comp**".**! with that 
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for AJ k FX GftB®, which was obtained previously by the same method J . The 
dots iu Fig. 2 represent the averaged time history of 48 APEX GRB» in the 
brood energy channel 100-1000 key. The two data sets agree quite well in spite 
of tha instrumental differences between APEX and BATsE. 

A comparison of “color* of GRB® along the averaged time history has been 
dour by separate averaging of their profiles for discriminator channels #2 aad 
if 3> aligned with respect to the brightest bins for the sum of those two channels. 
The average time profile of the hardness ratio ( HR), defined aa the ratio of counts 
m 100-300 k*V to the counts in 50-100 keV, is shown in Pig. 3. The HR ha* a 
maximum around 0.0 during the time interval from -0 to 0, during which the 
averaged time history rise® to its maximum from the lewd of 0.1. Before the rise 
front, i.c., within lime interval from -16 to -10, the averaged HR ia about 0.7. 
After the maximum, HR decreases with the averaged ftux to the value at the 
beginning of the rise front. Thus, we conclude that HR is u maximum during 
the rise front. 

Ill COMPARISON BETWEEN STRONG AND WEAK EVENTS 

The total set of GRB® mov be sub divided into “strong” and “weak” 
events. Wc first consider a simple separation, where strong wui weak events 
are associated with VfV mkX < 0.34 and > 0.34, respectively (for the sat of 
DATSE GRBs {V{V m9A ) = 0.34 ± 0.02 1 ). The averaged time historic® of 112 
strong and 79 weak event® thv» defined arc shown in Fig. 4 by bold and thin 
lines, respectively. No statistically significant difference ia seen between the two 
profiles. 

In cosmological model® the Urns dilation stretch factor between the stroug 
and weak event profile® should he ~ (1 + {j w ))/(1 + (x,}) t where {^w) and 
are the average redehift® of the weak and strong event®, respectively. Strictly 
speaking, the above comparison® arc compromised in a coamn logical model by 
the redshift of photon® from higher to lower energies. Because of the redshift, 
the profile of weak bursts in the observed energy channel corresponds to higher 
rest energies than the profile of strong bursts. Since the burst® axe intrinsically 
narrower at higher energies (Fig. 2), this effect could wipe out any time dilation 
which might be present. 

In or dot to compensate for this effect, we compared the averaged time 
history of 38 “very strong" event® (V7K*ui < 0.06) in 100-300 keV and the 
averaged time history of 42 “vary weak” events (Vfv mKM > 0.00) in 50-100 keV. 
Th« result® (Fig. 5) again show no evidence for systematic stretching of tlie 
very weak t*vauw icloiiv^ U 1AV tuj, y •imun t.n.utu F«r w n!«n 

®how in Fig. 5 the stretched profile expected from cosmological models whnift 
acceptably fit the BATSE munber/intenrity distribution 1 . Here wo assumed 
{V/Vmm*) ^ -026 for tho very strong events and (V/V«**) » 8 for the very wvak 
events, from which Fig. 2 of Ref. 7 predict* stretch factor® of 22-67% Fig. 5 
represents a conservative test since the rest energies uf the weak events are <x\ 
the average lower tliao those of the strong event® for these value® of (x), Thus, 
the redahift in energy of photon® cannot account for the absence of time dilation. 

On the other bund, the averaged HR tuna profile® show a difference l>e- 
twocn strong and weak subsets (Fig. 6). Both profiles have rut her similar shop** 
with a fiat maximum along the rising part of the averaged time histories and 
a dccrcoee just alter the brightest time bin, but the averaged values of HR for 
strong events nre larger than the corresponding averaged value" for weak events. 
A similar correlation Wwmi intenoity and hordn*** was found previously for 
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APEX GRBs 2 This result also confirms and extaods tha correlation found 
previously for BATSE GRBe*. 

TV. CONCLUSIONS 

1) The averaged time history of GRBe is non- symmetric: from the level of 
0.25 of its maximum height the cuie front lasts about 2—3 s and back slope 
about 4 *5 s. The averaged HR profile has a broad maximum along the 
rise front of the averaged time history, and decreases along its back slope. 

2) We find uo evidence in this analysis for time- dilation in the averaged time 
hii tones of ^weak** GRBs relative to ^vtcorng** avnmitK. Ff wri iiiiWUTno n. 
cosmological sotmwio, th« redshift of narrower temporal structure from 
higher to lower energies Is insufficient to account for the absence of a time 
dilation effect, 

3} We find evidence for differences between strong and weak CRBs in the 
averaged hardness ratio: strong events have larger HR then weak ones 
r»ot only averaged over time hut also at each point in tho averaged lime 
profile. In the absence of a ooantoiogical interpretation, we must then 
consider the similarities of the HR profiles to indicate their association 
with the same unique class of radiating objects, but with some important 
difference between them which leads to a hardness/ intensity correlation. 

FWther analysis is essential to compare the sensitivity of this method 
with that of other analyses which have been shown to support a coeinological 
interpretation 9 . We will conti n u e to investigate these issues further with various 
refinements to our analysis, including a larger sample of burets, vwu of data with 
fiinir time resolution, aud use of data in other energy ranges. 
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Pig, 1 '-Flux v«r»tia time. 




Fig. r> — Flux versus time. Bold line 
correspond to 41 very strong" events 
with innui end* stretched by 2 2% 
(vml outer eu«JU aticiched by 6T%. 
Thm line carru»|,KJiula to "very 
vfenk* events- 


Fi«. 2- Flux versus time. Bold liiue: 
50-100 loeV. Thin line: 100— CUOO 
keV. Dots represent APEX data. 




Fig. 0 — Hardness ratio versus time. 
Bold line correspond to “strong” 
ewmfcs. Thin Hne corresponds to 
events. 







